Observation of the L>sj(2317) and Dsj(2457) in B decays 
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We report the first observation of tlie B — > l)_Da,j(2317) and B — > DDsj (2457) decays based on 
123.8 X 10'^ BB events collected with the Belle detector at KEKB. We observe the Dsj{2317) decay 
to Dsit^ and -Dsj(2457) decay to the D*tt^ and Ds^ final states. We also set 90% CL upper limits for 
the decays _D«,7(2317) L>*7, L>,,j(2457) L>*7, D.j(2457) ~* D.7r° and D,j(2457) -» DsTt+tt". 

PACS numbers: 13.25.Hw, 14.40.Lb 



Recently a new DsTt'^ resonance with a mass of 
2317 MeV/c^ and a very narrow width was observed by 
the BaBar collaboration Q. A natural interpretation is 
that this is a P-wave cs quark state that is below the 
DK threshold, which accounts for the small width 
This interpretation is supported by the observation of a 
D*!:^ resonance ;3:| by the CLEO collaboration 1] and 
Belle collaboration ^] . All groups observe these states in 
inclusive e^e^ processes. The mass difference between 
the two observed states is consistent with the expected 
hyperfine splitting of the P-wave Dg meson doublet with 
total light-quark angular momentum j = 1/2 0. How- 
ever, the masses of these states are considerably below 
potential model expectations 0, and are nearly the same 
as those of the corresponding cu states recently measured 
by Belle 0. The low mass values have caused specula- 
tion that these states may be more exotic than a simple 
qq meson system H H El [H [H [3. To clarify the 
nature of these states, it is necessary to determine their 
quantum numbers and decay branching fractions, par- 
ticularly those for radiative decays. In this context it is 
useful to search for these states, which we refer to as Dsj, 
in exclusive B meson decay processes. 

We search for decays of the type B DDgj, which are 
expected to be the dominant exclusive Dsj production 
mechanism in B decays. Because of the known properties 
of the parent B meson, angular analyses of these decays 
can unambiguously determine the Dgj quantum num- 
bers. Moreover, since QCD sum rules in HQET predict 
that P-wave mesons with j = 1/2 should be more read- 
ily produced in B decays than mesons with j = 3/2 [31, 
the observation of P — > DD^j would provide additional 
support for the P-wave nature of these states as well as 



serving as a check of these predictions. 

In this Letter we report on a search for the B — > 
DDsj{2M7) and B DDsj{2A57) decays based on a 
sample of 123.8 x 10^ BB pairs produced at the KEKB 
asymmetric energy e~^e~ collider |l5|. The inclusion of 
charge conjugate states is implicit throughout this report. 

The Belle detector has been described elsewhere [l^ . 
Charged tracks are selected with a set of requirements 
based on the average hit residual and impact parameter 
relative to the interaction point (IP). The transverse mo- 
mentum of at least 0.05 GeV/c is required for each track 
in order to reduce the combinatorial background. 

For charged particle identification (PID), the com- 
bined information from specific ionization in the central 
drift chamber (dE/dx), time-of-flight scintillation coun- 
ters and aerogel Cerenkov counters is used. Charged 
kaons are selected with PID criteria that have an effi- 
ciency of 88%, a pion misidentification probability of 8%, 
and negligible contamination from protons. All charged 
tracks with PID responses consistent with a pion hypoth- 
esis that are not positively identified as electrons are con- 
sidered as pion candidates. 

Neutral kaons are reconstructed via the decay Kg — > 
TT+TT" with no PID requirements for the daughter pi- 
ous. The two-pion invariant mass is required to be within 
9 MeV/c^ (~ 3cr) of the K'^ mass and the displacement 
of the 7r"'"7r~ vertex from the IP in the transverse {r — Lp) 
plane is required to be between 0.2 cm and 20 cm. The 
direction in the r — Lp plane from the IP to the tt+tt^ 
vertex is required to agree within 0.2 radians with the 
combined momentum of the two pious. 

Photon candidates are selected from calorimeter show- 
ers not associated with charged tracks. An energy de- 
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position of at least 30 MeV and a photon-like shape are 
required for each candidate. A pair of photons with an 
invariant mass within 12 MeV/c^ (~ 2.5a) of the 7r° mass 
is considered as a candidate. 

We reconstruct D^{D^) mesons in the K^tt^ , 
K'^TT^TT^TT^ and K^tt^tt'^ (K^tt^tt^) decay channels 
and require the invariant mass to be within 12 MeV/c^ 
(1.5(7 for i<'+7r~7r'' and 2.5a for other modes) of the 
D°{D-) mass. For the 7r° from the D° K+tt'tt^ 
decay, we require that the tt" momentum in the T(45') 
center-of-mass (CM) frame be greater than 0.4 GeV/c in 
order to reduce combinatorial backgrounds. We recon- 
struct mesons in the (/'7r+, K*^K'^ and K^K^ decay 
channels. </) mesons are reconstructed from K'^K~ pairs 
with an invariant mass within 10 MeV/c^ (2.5r) of the 
(j) mass. K*^ mesons are reconstructed from K^tt^ pairs 
with an invariant mass within 75 MeV/c^ (l-5r) of the 
K*^ mass. After calculating the invariant mass of the 
corresponding set of particles, we define the D'^ signal re- 
gion as being within 12 MeV/c^ (~ 2.5a) of the Ds mass. 
D* mesons are reconstructed in the D* — > Ds^ decay 
channel. The mass difference between I?* and Dg candi- 
dates is required to be within 8 MeV/c^ of its nominal 
value (~ 2.5cr). The Dgj candidates are reconstructed 
from D^s"^ mesons and a tt", 7, or tt^tt" pair. The mass 
difference M{Dsj) - M{D^*^) is used to select J can- 
didates. We use central mass values of 2317 MeV/c^ and 
2460 MeV/c^ for Dsj{2?,n) and Dsj{2A57) respectively 
and define signal regions within 12 MeV/c^ for the cor- 
responding mass difference. 

We combine D and Dsj candidates to form B mesons. 
Candidate events are identified by their CM energy differ- 
ence, iS.E = (^^ Ei) — £^bcam, and the beam constrained 

mass, Afbc = V^bcam ~ {Y^iPiY^ where i^beam is the 
beam energy and pi and Ei are the momenta and ener- 
gies of the decay products of the B meson in the CM 
frame. We select events with 5.272 GeV/c^ < Mbc < 
5.288 GeV/c^ and |A£;| < 0.2 GeV, and define a B sig- 
nal region of |Ai?| < 0.03 GeV. In cases with more than 
one candidate in an event, the one with D and ' 
masses closest to the nominal values is chosen. We use a 
Monte Carlo (MC) simulation to model the response of 
the detector and determine the efficiency 

Variables that characterize the event topology are used 
to suppress background from the two-jet-like e^e~ —>■ qq 
continuum process. We require |cos0thi | < 0.80, where 
6'thr is the angle between the thrust axis of the B candi- 
date and that of the rest of the event; this eliminates 77% 
of the continuum background while retaining 78% of the 
signal events. To suppress combinatorial background we 
apply a restriction on the invariant mass of the D meson 
and the Tr° or 7 from Dsj decay: M{Dtt°) > 2.3 GeV/c^, 
M{D"f) > 2.2 GeV/c2. 

The AE and Dgj candidate's invariant mass (M(Ds,/)) 
distributions for B — s- DDgj candidates are presented 



in Fig. ^ where all -D° and D~ decay modes are com- 
bined. Each distribution is the projection of the sig- 
nal region of the other parameter; distributions for 
events in the M{Dgj) and AE sidebands are shown as 
crosshatched histograms. Clear signals are observed for 
the DD,j(2317)[n,7rO] and DD,j(2457)[D*7r°, D37] fi- 
nal states. The measured masses for the (2317) 
and Dsj{2A57) are (2319.8 ± 2.1 ± 2.0) MeV/c^ and 
(2459.2 ± 1.6 ± 2.0) MeV/c^ respectively. The fitted 
widths are consistent with those expected for Dsj mesons 
of zero intrinsic width. The systematic error in the Dgj 
mass is expected to come from the photon energy scale. 

We also study the helicity distribution for the 
Ds,/(2457) Dsj decay. The helicity angle Oos-y is de- 
fined as the angle between the -Dsj(2457) momentum in 
the B meson rest frame and the Ds momentum in the 
Dsj(2457) rest frame. The Oos-y distribution in the data 
(Fig- El is consistent with MC expectations for the J — 1 
hypothesis for the Z3sj(2457) (x^/n.d.f= 5/6), and con- 
tradicts the J = 2 hypothesis (x^/n.d.f.= 44/6). The 
J = hypothesis is already ruled out by the conservation 
of angular momentum and parity in Dgj {2A57) — > Dsj. 

For each decay channel, the AE distribution is fitted 
with a Gaussian signal and a linear background function. 
The Gaussian mean value and width are fixed to the val- 
ues from a MC simulation of signal events. The region 
AE < —0.07 GeV is excluded from the fit to avoid con- 
tributions from other B decays of the type B DDsjX 
where X denotes an additional particle that is not re- 
constructed. The M{Dsj) distribution is fitted by the 
sum of a Gaussian for the signal, and a linear func- 
tion for the background. The Gaussian width is fixed 
to the value found in the MC (6-7 MeV/c^ depending 
on the decay mode). The fit results are given in Tabled 
where the listed efficiencies include intermediate branch- 
ing fractions. We use the AE distribution to calculate 
the branching fractions. The statistical significance of the 
signal quoted in Table is defined as ^J —2h\{CQ / Cmax) , 
where Cmax and Cq denote the maximum likelihood with 
the nominal and with zero signal yield, respectively. 

The results of combined fits of B^ — > D^D'^j and 
B'^ D~ D~^j modes assuming isospin invariance are 
shown in Table UTI The normalization of the background 
in each sub-mode is allowed to fioat while the signal yields 
are required to satisfy the constraint Ni ~ N^g ■ B{B ^ 
DDsj) ■ Ei , where the branching fraction B{B DDgj) 
is a fit parameter; A^^^ is the number of BB pairs 
and Ei is the efficiency, which includes all intermediate 
branching fractions. From the two B DDsj (2457) 
branching fraction measurements, we determine the ra- 
tio BiDsji2A57) Ds-f)/B{Dsji2457) ^ /?>«) = 
0.38 ±0.11 ±0.04. 

_ The signals for the B DDsj{2317)[DsTr°] and B 
DDsj{24:57)[Dl'K^,Dsj] channels have greater than 5a 
statistical significance. Figure |31 shows the AE distribu- 
tions for the other channels, where significant signals are 
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FIG. 1: AE (left) and M(Dsj) (right) distributions for the B DDsj candidates: (a) Dsj(2317) ^ DsTt'^, (b) Dsj(2457) ^ 
D*7r° and (c) _Da.j(2457) ^ Ds-f. Points with errors represent the experimental data, crosshatched histograms show the 
sidebands and curves are the results of the fits. 



TABLE I: Product branching fractions for B — > DDsj decays. 



Decay channel 




AE yield 


M{Dsj) yield 


Efliciency, 10"* 


B, 10-* 


Significance 


B+ D°Df j{2317) 


[D+n% 


13.7lli 




1.36 


8.ll|^ ± 2.4 


5.0cr 


B° ^ D-D+^(2317) 


lDtn°] 


10.31^3:? 




0.97 


S.Gtli ± 2.6 


6.1a 


B+ D°D+j{2317) 


[Ds + -f] 


o 4+2.8 




1.08 


2.5ttl{< 7.6) 




B° D-Dfj(2S17) 




9 Q+2.5 


^■"-1.9 


0.69 


2.7l^:^(<9.5) 




B+ D°D+j{2A57) 




7 9+3.7 
' -^-3.0 


8.9tti 


0.49 


11.9llg ±3.6 


2.9a 


B° D-D+j{24:57) 




11.81^:1 


14.91*:^ 


0.42 


22.7le:2 ±6.8 


6.5a 


B+ D°Dfj{2'i57) 






20.21^:^ 


2.75 


5.6li;^ ± 1.7 


5.0a 


B° -* L>-D+j(2457) 


[Dt-/] 


18.511° 


19.611^ 


1.83 


8.2l^;2 ± 2.5 


6.5a 


B+ D°Dfj{2A57) 


[D's+y] 




8.2tt-l 


1.15 


3.ll^:^(<9.8) 




B° ~* L>-L>+j(2457) 




1 1+1.8 
^■^-1.2 


0.211:1 


0.71 


i.3l?:!^(< 6.0) 




B+ D°D+j{2457) 


[Dt7T+n-] 


< 4.0 


n + 2.0 
~^.^-1.6 


1.89 


< 2.2 




B" L>-D+^(2457) 


[L>+7r+7r-] 


< 2.5 


-1 2+2-^ 
r.^-2.0 


1.35 


< 2.0 




B+ D°Djj{2'i57) 




< 2.4 


i.ol^:^ 


0.94 


< 2.7 




B° ^ D- D+j{2457) 




< 2.4 


0.3li« 


0.68 


< 3.6 





not seen. We set 90% confidence level (CL) upper limits 
for these modes. 

We study the possible feed-across between all stud- 
ied Dsj decay modes using MC. We also analyze a MC 
sample of generic BB events corresponding to our data 
sample. No peaking background is found. As a check, 
we apply a similar procedure to decay chains with the 
similar final states: B Z)(*)_Di*'. For each mode, we 
measure branching fractions that are consistent with the 
world average values 



The following sources of systematic errors are consid- 
ered: tracking efficiency (1-2% per track), kaon identi- 
fication efficiency (1%), 7r° efficiency (6%), A'^ recon- 
struction efficiency (6%), D branching fraction uncer- 
tainties (2%-6%), signal and background shape param- 
eterization (4%) and MC statistics (3%). The uncer- 
tainty in the tracking efficiency is estimated using par- 
tially reconstructed D*~^ — > D^[K gTr~^TT~]Tr^ decays. 
The kaon identification uncertainty is determined from 
Z3*+ D'-^[K^n^]TT^ decays. The tt" reconstruction 
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TABLE II: Combined fit results. 



Decay channel 



Significance 



B -> DD, J {22,17) [D,7r°] 
B^DD, .,(2217) [Dl^\ 
B^DDsj(2457) [DI-k°] 
B DD, J (2457) 
B^L»D,j(2457) [Dl'y] 
B -> DD, J (2457) [DsTt+tt" 
B ^ DD,., (2457) [D^TfO] 



8.5t5;^ ± 2.6 
2.5t?:«(< 7.5) 

IT.Slg g ± 5.3 

6.7li;2 ± 2.0 
2.7t\%{< 7.3) 

< 1.6 

< 1.8 



6.1cr 
l.Scr 
6.4cr 
7.4cr 
2.1(7 




cos(eD3^) 

FIG. 2: The L>^,j(2457) ^ i3s7 helicity distribution. The 
points with error bars are the results of fits to the AE spectra 
for experimental events. Solid and dashed curves are MC 
predictions for the J = 1 and J = 2 hypotheses, respectively. 
The highest bin has no events because of the cut on the D'y 
invariant mass. 



uncertainty is obtained using decays to K~n'^ and 
K~'K^'K^ . We assume equal production rates for B^B^ 
and B^B° pairs and do not include the uncertainty re- 
lated to this assumption in the total systematic error. 
For the calculation of the branching fractions, the errors 
in the meson branching fractions are taken into ac- 
count. These uncertainties are dominated by the error on 
the Dg — > 4nT~^ branching ratio of 25% [11. The overall 
systematic uncertainty is 30%. 

In summary, we report the first observation of _B ^ 
Z)L»sj (2317) and B ^ DD^j (2A57) decays. The mea- 
sured branching fractions with the corresponding statis- 
tical significances are presented in Table UTI The observa- 
tion of Dsj(2457) — > Dsj decay eliminates the zero spin 
of Dsj(2457). The angular analysis of this decay sup- 
ports the hypothesis that the Dsj{2457) is a 1+ state. 

We wish to thank the KEKB accelerator group for the 
excellent operation of the KEKB accelerator. We ac- 




FIG. 3; AE distributions for the decay channels with not sig- 
nificant signals: (a) _DD,j(2317)[D*7], (b) Di3,,7(2457)[L>*7], 
(c) L>D,j(2457)[7:»,7r+7r-], (d) 73L»^j(2457)[D,7r°]. Open his- 
tograms represent the experimental data and curves show the 
results of the fits. 
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